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ABSTRACT 

In this paper, the effect of the intrinsic distribution of cosmological candles is investigated. 
We find that, in the case of a narrow distribution, the deviation of the observed modulus of 
sources from the expected central value could be estimated within a ceratin range. We thus 
introduce a lower and upper limits of x 2 , Xmin anc ^ Xmax' t0 estimate cosmological parameters 
by applying the conventional minimizing x method. We apply this method to a gamma-ray 
burst (GRB) sample as well as to a combined sample including this GRB sample and an SN la 
sample. Our analysis shows that: a) in the case of assuming an intrinsic distribution of candles 
of the GRB sample, the effect of the distribution is obvious and should not be neglected; b) 
taking into account this effect would lead to a poorer constraint of the cosmological parame- 
ter ranges. The analysis suggests that in the attempt of constraining the cosmological model 
with current GRB samples, the results tend to be worse than what previously thought if the 
mentioned intrinsic distribution does exist. 
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1 INTRODUCTION 

One of the greatest achievements obtained in the past few years in 
astrophysics is the determination of cosmological parameters with 
type la supernovae (SN la), which suggests an accelerating universe 
at large scales (Riess et al. 1998, Perlmutter et al. 1999, Tonry et 
al. 2003, Barris et al. 2004, Knop et al. 2003, Riess et al. 2004). 
The cosmic acceleration was also confirmed, independently of the 
SN la magnitude-redshift relation, by the observations of the cos- 
mic microwave background anisotropies (WMAP: Bennett et al. 
2003) and the large scale structure in the distribution of galaxies 
(SDSS: Tegmark et al. 2004a, 2004b). It is well known that all 
known types of matter with positive pressure generate attractive 
forces and decelerate the expansion of the universe. Given this, a 
dark energy component with negative pressure was generally sug- 
gested to be the invisible fuel that drives the current acceleration 
of the universe. There are a huge number of candidates for the dark 
energy component in the literature, such as a cosmological constant 
A (Carroll et al. 1992), an evolving scalar field (referred to by some 
as quintessence: Ratra and Peebles 1988; Caldwell et al. 1998), the 
phantom energy, in which the sum of the pressure and energy den- 
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sity is negative (Caldwell 2002), the so-called "X-matter" (Turner 
and White 1997; Zhu 1998; Zhu, Fujimoto and Tatsumi 2001; Zhu, 
Fujimoto and He 2004b), the Chaplygin gas (Kamenshchik et al. 
2001; Bento et al. 2002; Zhu 2004), the Cardassion model (Freese 
and Lewis 2002; Zhu and Fujimoto 2002, 2003, 2004; Zhu, Fu- 
jimoto and He 2004a), and the brane world model (Randall and 
Sundrum 1999a, 1999b; Deffayet, Dvali and Gabadadze 2002). 

Samples of SN la sources available in the early analysis con- 
tain only sources with redshifts z < 1. Although observations 
of the fluctuation in the cosmic microwave background (CMB) 
can constrain the cosmological model up to redshifts as high as 
2 ~ 1000 (e.g., Spergel et al. 2003), a more direct measurement of 
the universe with objects located at very large distances is strongly 
desired. Fortunately, recent observations extended the SN la sam- 
ple to sources with redshifts as large as z = 1.7. The previous 
result was confirmed by these high redshift sources and the anal- 
ysis revealed that before its acceleration the universe underwent a 
period of deceleration (Riess et al. 2004). The success of includ- 
ing high redshift SN la sources inspires us to great efforts to search 
for cosmological rulers with much higher redshifts. Based on the 
E p — E~ t relation found recently in a class of gamma-ray bursts 
(GRBs) (Ghirlanda et al. 2004b), Dai et al. (2004) assumed that the 
GRB sources obeying this relation can be used to measure the uni- 



2 Qin et al 



verse. In their sample of 12 GRBs, two have redshifts z > 2. Soon 
after their work, the same issue was investigated by many authors 
(see Ghirlanda et al. 2004a; Friedman and Bloom 2005; Firmani 
et al. 2005; Xu et al. 2005; Liang and Zhang 2005). It was found 
that current GRB data which are lack of low redshift sources could 
be used to marginalize some parameters in their reasonable ranges 
(see Xu et al. 2005 for a detailed explanation), or they could be em- 
ployed to constrain the cosmological model with a new Bayesian 
method (Firmani et al. 2005). Although the size of the current GRB 
sample is small and low reshift sources are missed, the idea that 
some high redshift extragalactic sources other than SN la might be 
employed to determine the cosmological model is quite interesting 
and promising. 

It would be natural that, for a kind of source which could 
serve as candles, one assumes a distribution of luminosity, which 
is reasonable due to fluctuation. As discussed in Kim et al. (2004), 
the uncertainty of a source must include both the systematic uncer- 
tainty and the magnitude dispersion. We argue that, if there exists 
a distribution of luminosity of the candles, the expected luminosity 
itself (or the corresponding deduced luminosity distance) could be 
different from source to source, which would be due to an intrin- 
sic property rather than to the measurement uncertainty. This raises 
a topic of finding an appropriate method to estimate cosmological 
parameter ranges with candles with a certain distribution. 

When employing candles such as SN la or GRBs to measure 
the universe, the confidence level associated with the fit of the the- 
oretical curve to the luminosity distance data was described by a 
statistic \ 2 which is defined under the assumption that the mea- 
surement uncertainty is the only cause of the deviation of the data 
to the curve. The best fit will be obtained when one reaches the 
minimum value of % 2 ■ However, for candles with a certain distri- 
bution, the deviation of the observed luminosity from the expected 
curve must be caused by both the measurement uncertainty and the 
distribution itself. When taking into account the distribution of lu- 
minosity, the x 2 statistic could not be defined if the distribution it- 
self is unknown. The minimizing \ 2 method will not be applicable 
if the statistic itself cannot be defined. 

In the following, we will study how to deal with this matter 
and investigate what one can expect from the analysis. A corre- 
sponding method will be proposed and will be illustrated with two 
samples. 



2 THE METHOD 

In this section, we propose a method to deal with candles with a 
certain distribution when employing them to constrain the cosmo- 
logical model. As mentioned above, the statistic \ 2 could not be 
defined for candles with a certain distribution if the distribution it- 
self is unknown. Even if the distribution is known, the statistic is 
still undefinable since there is no way to know the real luminos- 
ity of each source. These difficulties lead to two problems. One is 
that the well-known minimizing \ 2 method could not be applicable 
without a definition of the statistic. The other is that the probabil- 
ity associated with the statistic x 2 , if we define it when taking into 
account the deviation arising from the distribution, is not available 
(since the real luminosity of each source is unknown). 

It is known that the convolution of two Gaussian is still a 
Gaussian with a width that is given by the quadratic sum of the 
two widths of the original distributions. That is a 2 = a\ + a\, 
where a\ and a 2 are the variances of the two Gaussian functions 
concerned and a 2 is that of the resulted Gaussian. 



Let us consider the deviation of an observed luminosity dis- 
tance modulus, fi b, of a source from the real value of the quantity, 
Hth, which follows 



(fiob ± <Job) - tJ-th(z; Hp, Q m , Qa) 
= (p b ± o ob ) - [fith,o(z; Hp, Q m ,£l A ) + A-Hth], 



(1) 



where a ob is the measurement uncertainty of /j, ob , ^th,p is the cen- 
tral value of fi th , which is the real value of the modulus expected 
in the case when there is no distribution of the candles, and A/xth 
represents the deviation of n t h from fith,o- Suppose that the distri- 
bution of candles is narrow enough so that the absolute value of the 
deviation of nth from fith,p, \A^, t h\, is small. According to the er- 
ror transform formula, the uncertainty of fi ob relative to fi thi p could 
be determined by 



>ob,p — 



(2) 



Relative to the expected central moduli, the x 2 statistic of a sample 
of the candles could be determined by 

[fJ.ob,i — fJ-th,o,i(z; Hp, Q m , &a)] 2 



x 2 = E ] 



(3) 



[Note that, in the case of SN la, a 2 ob i should be replaced by 
a ob,i + a vi where a b,i is the uncertainty in the individual distance 
moduli deduced from the empirical relation between the light-curve 
shape and luminosity and a v is the uncertainty associated with the 
dispersion in supernovae redshift (transformed to units of distance 
moduli) due to peculiar velocities (see Riess et al. 2004)] 

It seems that, with equation (3), one might be able to evalu- 
ate the x 2 statistic. But because Afi thii is in no way to be known, 
this is unfortunately not true. However, under the condition that the 
distribution of candles is narrow, we can estimate Afi t h.i with the 
width of the distribution. Let ajis be the width of the distribution of 
lithj ' iHhfi (called the intrinsic distribution of the relative luminos- 
ity distance moduli). (Note that fi th /fi t h,o should of course become 
unity when there is no deviation of n t h from fith.o)- We assume 
\Afi t h,i\ — (Tdis^thfi,i- Thus the x 2 statistic could be estimated by 



X 



E 



[^ob,i — fJ-th,o,i{z; Hp, Q m , Qa)] 2 



'ob, 



(4) 



As long as a dis is provided, the x' 2 statistic is then avail- 
able according to (4). For any kind of candle, quantity Odis could 
be estimated when the sample employed is large enough and the 
measurement uncertainty a ob is small enough and when the cos- 
mological model is fixed. Obviously, this could not be realized at 
present since the cosmological model itself is currently a target to 
be pursued and for interesting candles the measurement uncertainty 
is always quite large. But this cannot prevent one to estimate the 
limits of adis- As the deviation of /x D 6 from Hth,o is caused by 
both the distribution of \i t h and the measurement uncertainty of 
LL ob itself, a dis must be smaller than <7 diSimax , where (j diSjmax is 
the width of the distribution of Hob/ l-i-th,p, which is determined 

by Sdi^max = n2{^ob,i/ W,o,i - I) 2 / (N - 1). with N be- 



ing the size of the sample. Let us over estimate the effect of the 
measurement uncertainty in the opposite way. Within the range of 
[/ioM — <Job,i, Hob,i + foi),i] we take the value that is the closest one 
to Hth,o,i as n* ob>i . Obviously, the distribution of [i ob j 'fi t h,o would 
be narrower than the distribution of fith/^th.o since the deviation 
caused by the measurement uncertainty is over subtracted. We take 
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the width of the distribution of y,* ob / ' Hthfi as a diStrnin , which is 
calculated with a d i S , min = 



E(m* m Mmm- i)7(^-i)- 



Clearly, Odis must be larger than With these two quan- 

tities we have 



and 



2 

Xmin 



2 

Xmax 



E 



/j2 _i_ ™-2 j/2 

ob.i 1 dis,maxr"£h,0,i 



^ob,i dis,minf-^th,0,i 



(6) 



Since (7d; s ,min < (7 d ; s < (7d; s ,max, One getS Xmin < X 2 < Xmax- 

With equations (5) and (6), one can calculate the corresponding 
probability associated with the x 2 statistic and confine the con- 
ventional confidence contour. In this way, cosmological parameters 
would be constrained. With this estimating method, the first prob- 
lem is largely eased and the second is solved. 



3 APPLICATION 

Let us consider a GRB sample. The sample was presented and stud- 
ied in Xu et al. (2005) and Xu (2005) (the XDL GRB sample) which 
contains 17 GRBs. As suggested in Ghirlanda et al. (2004a), the 
scatter of the data points of their GRB sample around the correla- 
tion of E p — found recently (Ghirlanda et al. 2004b) is of a very 
small order. 

To check if the data of the XDL GRB sample are consis- 
tent with no scatter beyond the measurement errors in terms of 
statistics, the simplest method is to calculate the mean of the de- 
viation of the deduced luminosity distance moduli from the ex- 
pected one of the sample and then compare it with the average 
of the measu rement error. The mean of the deviat ion is defined 
as a dev = J2((^ob,i - Hex,i)/Hex,i) 2 /(N - 1), where fi ex is 



the expected value of /j,, while the average of the measurement 
error is calculated with a err = ^Jj2( a ob,i/Hex,i) 2 /(N - 1). 

(Note that, as redshifts of these sources are not the same, we con- 
sider the relative values.) We get the following from the XDL 
sample: adev = 0.0122 and a err = 0.0116, where we adopt 
(fi m , fi A , h) = (0.29, 0.71, 0.65). It shows that the deviation is 
slightly larger than the measurement error. (Ignoring the slight dif- 
ference between the two quantities, the result confirms what sug- 
gested in Ghirlanda et al. 2004a, 2004b.) Taking Hthfi as n ex 
adopted here, one finds that a dev is identical with <7d» s>max de- 
fined in last section. Thus, for the XDL sample, adis < 0.0122, 
suggesting that the distribution, if exists, would be quite narrow. 
Another approach involves a simulation analysis. We assume that 
there is no intrinsic distribution of the deduced luminosity distance 
moduli, and thus the deviation observed is due to the measurement 
uncertainty. Obviously, under this assumption the distribution of 
Hob/ Hex should peak at unity. According to the null hypothesis, 
the observed value of Hex for each source is obtained by chance 
from a parent population of Hob whose distribution obeys a Gaus- 
sian with the measurement uncertainty served as the width of the 
Gaussian. For each source one can create a y! ob via simulation as 
long as the expected value Hex and the measurement uncertainty 
are known. In this way, from the 17 Hex and the corresponding 
measurement uncertainties, one can create a set of 17 Hob data 
by a Monte-Carlo simulation and then obtain a set of 17 Hob/ Hex 



data. We perform 100 times of simulation and get 100 sets of 17 
Hob/ Hex data. Combining these 100 sets we get a large sample with 
its size being 1700. The deviation of the relative simulated luminos- 
ity distance moduli from the expected one (the unity) is defined as 

v'dev = ( /EW M Ki-l) 2 /(JV-l). Note that a' dev could 



be written as a' dev 



Y.((H' obti -Hex,i)/Hex,i) 2 /(N-l), 



which could thus be directly compared with adev, the deviation 
of the observed data defined above. From the XDL sample we get 
a' dev — 0.0113, which suggests that the deviation associated with 
observation, denoted by a dev , is also slightly larger than that ex- 
pected from the measurement uncertainties. Two methods come to 
almost the same result, suggesting that there might be an intrinsic 
distribution of the relative luminosity distance moduli of the XDL 
sample, although it would be quite narrow (as the difference be- 
tween a'de V an d o- dev and that between a err and a dev are small). 

To illustrate how to apply the method proposed above to deal 
with data with intrinsic distributions, we assume in the following 
that there is a distribution of the true value of the deduced rela- 
tive luminosity distance moduli for the XDL sample, although the 
distribution, if it exists, might be very narrow (see what suggested 
above). For the sake of comparison, we perform the fit with three 
X 2 statistics. One is the conventional x 2 which could be determined 
by (3) when taking Afi t h,i = 0. The other two are Xmin and Xmax 
which are determined by equations (5) and (6) respectively. Each 
X 2 statistic is calculated with the XDL GRB sample in many tries. 
In each try, we adopt a set of parameters and based on these pa- 
rameters we deduce both the observed and theoretical luminosity 
distance moduli. With these moduli and the measurement uncer- 
tainties, we are able to evaluate a d i S ,min and (Tdis^ax (see what 
proposed in last section), and then the corresponding x 2 statistic 
would be well determined (Ho = 65kms~ 1 Mpc~ x is adopted 
throughout this paper). For each x 2 > the best fit will be obtained 
when the smallest value is reached. 

Displayed in Fig. 1 are the Hubble diagram and the confident 
contour plot of the XDL GRB sample. As concluded previously 
by other authors (see Ghirlanda et al. 2004a; Friedman and Bloom 
2005; Xu et al. 2005), currently, employing GRB samples alone 
cannot tightly constrain the cosmological model. Fig. 1 shows that, 
the parameter ranges are indeed poorly constrained even there is no 
intrinsic distribution of the relative luminosity distance moduli (see 
solid lines in Fig. lb). Taking into account an intrinsic distribution 
of the moduli leads to much poorer results. This indicates that if 
there indeed exists an intrinsic distribution of the moduli, the effect 
arising from the distribution should not be ignored. 

Shown in Table 1 are the best fit cosmological parameters for 
the three kinds of universe, obtained by applying the minimizing x 2 
method to the three x 2 statistics, where the la errors are estimated 
from the corresponding la contours in Fig. lb. As shown in Fig. 
lb, the la contours are not closed within the ranges of the plot. 
This leads a poor constraint to the limits of the parameters. Some 
limits are therefore not able to be determined, which are denoted 
by"?" in Table 1. 

The fact that the parameter ranges are poorly constrained 
(even when the intrinsic distribution of the relative luminosity dis- 
tance moduli is ignored) might probably be due to the lack of low 
redhsift sources, as it is already known that low redhsift sources are 
important when employing a GRB sample to constrain the cosmo- 
logical parameters (see Firmani et al. 2005). We thus follow what 
were done previously (see Ghirlanda et al. 2004a) to combine an 
SN la sample and the XDL sample to constrain the cosmological 
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Table 1. Best-fit parameters obtained by the least square method 



Sample 


Universe 


(n M ,n A , x §,„)° 




( n M . a A 


X 2 ■ ) 


SN+GRB 


flat 


(0.283+5J'q288 - 0.717, 197.9) 


(0.288 + q q201 , 0.712, 193.1) 


(0.288+^ 


SoolVs- - 712 ' 185 - 4 ) 


SN+GRB 


open 


(0.368 j:°^,0.857+°;^J, 197.0) 


,„,„, +0.0201 +0.0198 . 
<°- 281 -0.0201 • °- 717 -0.0353 ' 193 -2) 


(0.28ll° 


00669 n 7 , 7 + 0.0100 . 
0134 • °- 717 -0.0201 1 185 - 4 ) 


SN+GRB 


closed 


/„ 9R1 +0.0334 n _, -+0.0296 ,„, m 
(0.281_ 0469 , 0.717_ Q 0gf)4 , 198.0) 


(0.428 + 01 47 ,0.942+°- 22 |, 191.3) 


(0.441 + 2 


147 n 967+0- 226 isi l 1 
0201 .0 967_ 03S1 , 183.1) 


GRB 


flat 


(0.188+°' j™, 0.812, 19.69) 


(0.188+ ? - 176 , 0.812, 15.14) 


(0.188+^- 


539 , 0.812, 7.43) 


GRB 


open 


(0.187 + 9/ 221 , 0.682 + "- 221 , 19.67) 


(0.187+^ , 0.256 + ^ , 14.96) 


+ ? 

(0.154^ 


, 0.391 + ^ , 7.40) 


GRB 


closed 


(0.187lS;??i. 0.817+2;^, 19.67, 


(0.187+ J/ 177 , 0.817+9/ 551 , 15.14) 


(0.187+ 1/ 


54 ,0.817±°/ 59 ,7.44) 



Xq „ is the reduced \ 2 calculated with equation (4) when assigning a^is = 0. 



model. The SN la sample employed is that presented in Riess et al. 
(2004) (the so-called gold set of SN la) which contains 157 sources 
(where, many low redshifts sources are included). In the same way 
and for the same reason we apply the minimizing x 2 method to 
the three x 2 statistics to find the best fit cosmological parameters. 
Note that, unlike what is shown in the case of the GRB sample, the 
deduced luminosity distance moduli of the SN la sources do not 
depend on the adopted cosmological parameters. 

It is known that, in estimating the deduced luminosity distance 
moduli of the SN la sources, deviations caused by different magni- 
tudes of the peak luminosity of the sources have been checked. In- 
deed, we find that the distribution of the relative luminosity distance 
moduli of the SN la sample is very narrow (the figure is omitted). 
This suggests that, if it still exists (possibly caused by the small de- 
viation from the adopted empirical relation between the light-curve 
shape and luminosity), the intrinsic distribution must be extremely 
narrow. Thus we ignore the intrinsic distribution of the relative lu- 
minosity distance moduli of the SN la sample, and consider only 
&dis, min and ^dis,max for the GRB sample when we calculate the 
corresponding Xmin an d Xmax for the combined sample (including 
the XDL GRB sample and the gold SN la sample). 

Shown in Fig. 2 are the Hubble diagram and the confident con- 
tour plot of the combined sample. One finds that, including the SN 
la sample significantly improves the constraint of the ranges of cos- 
mological parameters. Once more, the result shows that taking into 
account the intrinsic distribution of the relative luminosity distance 
moduli leads to a poorer constraint. The effect is still obvious (al- 
though it is less obvious than that adopting the GRB sample alone) 
and therefore should not be neglected. 

The resulting best fit cosmological parameters as well as their 
lcr errors are listed in Table 1 as well. 



4 DISCUSSION AND CONCLUSIONS 

The effect of the intrinsic distribution of cosmological candles is 
investigated in this paper. Due to fluctuation, it is natural that a 
property (say, the luminosity) of sources served as a cosmological 
candle might form a distribution and scatter around a central value. 
If the distribution does exist, the statistic x 2 cannot be defined since 
the distribution itself is unclear and the real value of the property 
for each source is unknown. However, when the distribution is nar- 
row, the deviation of the observed modulus of each source from 
the central value could be estimated within a ceratin range. We ac- 



cordingly define a lower and upper limits of x^ 



' Xmin 



and Xn 



, to 



estimate cosmological parameters via the conventional minimizing 
X 2 method. The confidence contours of these two x 2 statistics can 
then be plotted in the conventional way, and with these curves the 
ranges of the parameters could be determined as long as a confi- 
dence level is assigned. 

With this method, a sample bearing a relatively small width of 
the intrinsic distribution of the deduced relative luminosity distance 
moduli would be applicable to constraining the cosmological pa- 
rameters. To illustrate this method we employ a GRB sample alone 
and later combine this GRB sample with the gold SN la sample, 
assuming that this GRB sample (the XDL sample) has an intrin- 
sic distribution of the deduced relative luminosity distance moduli 
while the SN la sample has not. The analysis suggests that: a) the 
effect of the intrinsic distribution of the relative luminosity distance 
moduli is obvious and therefore should not be neglected if the dis- 
tribution itself does exist; b) taking into account this effect would 
lead to a poorer constraint of the ranges of cosmological parame- 
ters. This indicates that in the attempt of constraining the cosmo- 
logical model with GRB samples, the results tend to be worse than 
what previously thought if the mentioned intrinsic distribution ex- 
ists, although the distribution is very narrow. 

As revealed recently by Wang et al. (2005), there is a clear evi- 
dence for a tight linear correlation between peak luminosities of SN 
la and their B — V colors at ~ 12 days after the B maximum. They 
found that this empirical correlation allows one to reduce scatters in 
estimating their peak luminosities from ~ 0.5 mag to the levels of 
0.18 and 0.12 mag in the V and I bands, respectively. We wonder 
if taking into account this effect can reduce the measurement un- 
certainty of the luminosity distance of the SN la sources. If so, the 
ranges of the cosmological parameters might be better constrained 
(when compared with Fig. 2) (this will be investigated later). 

As encountered in other cases, our method suffers from possi- 
ble evolution of candles. Quite recently, Firmani et al. (2004) found 
evidence supporting an evolving luminosity function of long GRBs, 
where the luminosity scales as (1 + z) 1 2 . It is unclear if the 
corrected gamma-ray energy, from which the luminosity distance 
moduli of the adopted GRB sample are deduced, evolves with red- 
shif. If so, the question if the GRB sample can still be used to con- 
strain the cosmological model should be answered. This deserves a 
detailed investigation. (It could be done only when the size of the 
sample is large enough). 
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Figure 1. (a) Hubble diagram of the XDL GRB sample, where the two edge 
curves (the dotted and dashed lines) of the flat universe are also represented, 
(b) Confident contour plot for the XDL GRB sample, where three dashed 
lines from the innermost curves to the outmost one represent the 1,2, and 3 
a levels of confidence calculated with the statistic Xmin respectively, while 
the three dotted lines represent those associated with the statistic Xmax 
respectively. For the sake of comparison, the confidence levels calculated 
without considering the distribution of Ji b/ Hth are also plotted (the solid 
lines). The straight line denotes the flat universe and the plus represents the 
best fit parameters of the flat universe obtained by the conventional mini- 
mizing x 2 method. 



Figure 2. (a) Hubble diagram of the combined sample, where the empty 
circle represents the XDL GRB sample and the filled circle stands for the 
gold SN la sample. Symbols of lines are the same as those denoted in panel 
(a) of Fig. 1 . (b) Confident contour plot for the combined sample, where the 
symbols are the same as those denoted in panel (b) of Fig. 1 . 
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